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Abstract. We derive new stability criteria for purely 
MHD instabilities in rotating jets, in the framework of the 
ballooning ordering expansion. Quite unexpectedly, they 
involve a term which is linear in the magnetic shear. This 
implies that cylindrical configurations can be destabilized 
by a negative magnetic shear as well as by a favorable 
equilibrium pressure gradient, in distinction with the pre- 
dictions of Suydam's stability criterion, which suggests on 
the contrary that the shear is always stabilizing. 

We have used these criteria to establish sufficient con- 
ditions for instability. In particular, the magnetic shear 
can always destabilize jets with vanishing current density 
on the axis, a feature which is generically found in jets 
which are launched from an accretion disk. We also show 
that standard nonrotating jet models (where the toroidal 
field dominates the poloidal one), which are known to be 
unstable, are not stabilized by rotation, unless the plasma 
(3 parameter and the strength of the rotation forces are 
both close to the limit allowed by the condition of radial 
equilibrium. 

The new magnetic shear-driven instability found in 
this paper, as well as the more conventional pressure- 
driven instability, might provide us with a potential energy 
source for the particle acceleration mechanisms underlying 
the high energy emission which takes place in the interior 
ofAGNjets. 

Key words: Magnetohydrodynamics : stability, balloon- 
ing, interchange - Jets : stability 



1. Introduction 

Jets are commonly observed in connection with Ac- 
tive Galactic Nuclei (AGNs) and Young Stellar Objects 
(YSOs). These jets are cylindrically collimated over re- 
markably long distances in comparison with their radial 
extents; the confinement provided by the magnetic ten- 
sion of the toroidal field in current carrying jets has long 
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been recognized as an efficient way to achieve this type of 
collimation (as shown, e.g., by Chan and Henriksen 1981, 
Blandford and Payne 1982, after an initial suggestion by 
Bcnford 1978). However, it is well-known from the field 
of thermonuclear fusion that toroidally confined plasma 
configurations are generically unstable, and instabilities 
in this context are able to destroy these configurations on 
very short time-scales (for a general background on MHD 
instabilities, in particular in the context of thermonuclear 
fusion, see Bateman 1978, and Freidberg 1987). Further- 
more, jets are subject to other instabilities, most notably 
the Kelvin-Hclmholtz one, which pose similar threats to 
their survival. On the other hand, there is ample evidence 
of instability in the observed jets, both directly (lateral 
displacements of the jet beam, bright knots...) and indi- 
rectly (models of jet synchrotron emission rely on particle 
acceleration through e.g. shock waves or turbulence), but 
the observed jet survival implies that (for reasons which 
are still unclear), instabilities in real jets lead to internal 
reorganisation and turbulence rather than to disruption. 

Most of the literature on jet stability has focused 
on the Kelvin-Hclmholtz instability (see, e.g., Birkinshaw 
1991 and references therein), both for hydrodynamical and 
MHD jets, and for a variety of jet equilibria. Two types 
of modes are produced by this instability: ordinary sur- 
face modes and reflected body modes. Surface waves are 
effectively confined to the jet interface with the external 
medium, while body modes appear and become dominant 
only for beam velocities in excess of the fast magnetosonic 
velocity (up to a factor of order unity). Modes growth rates 
decrease when the Mach number increases. A longitudinal 
magnetic field has a stabilizing influence due to magnetic 
tension; in particular, sub-Alfvenic (up to a factor of or- 
der unity) flows are completely stabilized. Radiative ef- 
fects can either enhance or reduce the growth rate of the 
instability depending on the steepness of the temperature 
dependence of the cooling function; in the presence of ra- 
diative cooling, the surface waves are apparently the most 
dangerous for jet survival (Hardee and Stone 1997; Stone 
el al. 1997). 
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Besides the Kclvin-Helmholtz instability, driven by 
velocity gradients, jets can be unstable with respect to 
purely MHD processes. In the context of ideal MHD, these 
instabilities arc usually divided into pressure-driven and 
current-driven instabilities (see e.g. Freidberg 1987 for de- 
tails). Pressure-driven instabilities are related to gradi- 
ents of the equilibrium pressure and to magnetic field 
line curvature. The excited modes include the so-called 
saussage mode, and are usually divided into interchange 
modes and ballooning modes; they share some common 
features with the Rayleigh- Taylor and Parker instabilities. 
Current-driven modes originate in the current parallel to 
the magnetic field, and include the so-called kink instabil- 
ities; the most dangerous is usually the m = 1 kink mode 
(m being the azimuthal wave- number) . 

Much less attention has been devoted to the analy- 
sis of these MHD instabilities, most probably because in 
superfast jets (whose beam velocity exceeds the fast mag- 
netosonic one), the kinetic energy exceeds the magnetic 
energy, and the Kelvin-Helmholtz instability is expected 
to dominate. Indeed, this expectation is borne out both in 
superfast and transmagnetosonic jets (superalfvcnic, but 
subfast beam velocities) for current-driven instabilities, 
whose growth rates are always substantially smaller than 
for the Kelvin-Helmholtz instability, at least for force-free 
jet equilibrium configurations (Appl and Camenzind 1992; 
Appl 1996). These authors also show that the jet current 
is stabilizing for supermagnetosonic flows, whereas it is 
destabilizing for transmagnetosonic ones. 

This paper focuses on pressure- and magnetic shear- 
driven instabilities, in the framework of the ballooning or- 
dering expansion. Pressure instabilities have been most ac- 
tively studied in fusion research (Kadomtsev 1966; Coppi 
et al. 1979; Dewar and Glasser 1983; see Freidberg 1987, 
and references therein). In astrophysics, they have been 
considered in the context of solar (e.g., Hood, 1986) and 
magnctospheric physics (see Ferriere et al. 1999 for a re- 
cent and synthetic overview of the situation in this field) . 
Although these instabilities have been virtually ignored 
in the context of MHD jets (with the exception of Begel- 
man 1998), it is important to characterize their conditions 
onset there; indeed, in spite of their local nature, they 
are able to produce large-scale disruptions over a dynami- 
cal time-scale in fusion devices. Furthermore, our analysis 
is also motivated by the following considerations. Small- 
scale Kelvin-Helmholtz and current-driven instabilities are 
stabilized by magnetic tension. Furthermore, the Kelvin- 
Helmholtz instability in the transmagnetosonic regime is 
confined to the jet boundary. These features don't make 
them favorable energy sources for the high energy emis- 
sion which takes place in jet interiors. By contrast, the 
pressure- and magnetic shear-driven instabilities consid- 
ered here usually have a small-scale component transverse 
to the magnetic field lines, and are expected to occur in 
the type of magnetic configurations which are typically 
considered for jet inner regions. 



In this paper, we clarify somewhat the relation between 
the various classical stability criteria found in the litera- 
ture, and we assess the role of the magnetic shear and of 
the rotation of the jet on the onset of jet instability, in 
particular in jet interiors. In the process, we point out the 
possible and unexpected destabilizing role of the magnetic 
shear. This paper is organized as follows. In Section 2, we 
recall the MHD equations and introduce our notations; 
we also briefly recall the origin of the ballooning ordering 
expansion, and present the perturbation equations in this 
approximation (their derivation is performed in the Ap- 
pendix); finally, we introduce the type of modes we analyse 
in this paper (which are interchange and not ballooning 
modes) and specialize the ballooning equations to cylindri- 
cal systems. In Section 3, we focus on the mode dispersion 
relation and discuss its properties for non-rotating and ro- 
tating jets. Section 4 discusses our results and concludes 
this paper. 

2. MHD equations and ballooning ordering 

For simplicity, we consider cylindrically symmetric jets, 
and we neglect the shear of both the vertical and angular 
velocities. We do not consider the question of the radial 
jet boundary. Quantities pertaining to the jet equilibrium 
configuration are labelled with a "0" subscript. 

2.1. MHD equations for jet equilibrium and perturbation 

In a frame which is rotating and moving with the jet, the 
MHD momentum equation reads 



dv VP, BVB , 

— = -v-Vv 1 h fl re - 2fl x v, 

at p p (i ' 

(1) 

where v is the fluid velocity in the moving frame, P* and p 
the fluid total (gaseous P and magnetic P m ) pressure and 
density, B the magnetic field, fi the vacuum permeabil- 
ity, r the cylindrical radius and e r the radial unit vector. 
This equation is closed as usual with the mass continu- 
ity equation, the induction equation in the flux freezing 
approximation, and an adiabatic equation of state (we as- 
sume that Pp" 1 is constant throughout the plasma for 
simplicity). 

The cylindrical equilibria considered in this paper are 
best characterized by introducing a number of vectors, 
K c = e /r Ve /h K b = VB /B , and K p = Wp /p , where 
e// = B Q / B Q is the unit vector in the direction of the un- 
perturbed magnetic field B ; K c is the curvature vector 
of the magnetic field lines, and Kb characterizes the in- 
verse of the spatial scale of variation of the magnetic field, 
while K p characterizes the inverse scale of variation of the 
fluid density. We also introduce the plasma (5 parameter: 
f3 = (Cs/Va) 2 where C$ and Va are the sound and Alfvcn 
speed. This parameter measures the relative importance 
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of the gas and magnetic pressures (our definition differs 
from the standard one by a factor 7/2). With these defini- 
titions, the jet force equilibrium relation reads 



PoVlflKp = p V\ (K c - 7C 6j _) + p Vt 2 re r . 



(2) 



This equation implies in particular that the component of 
K p parallel to the magnetic field vanishes. 

Linearizing the momentum equation Eq. (pi) , we obtain 



ddo 
~dt 



V(5P» (B -V)6B {5B-V)B 



(h 



PoP 



5 P n 2 re r - 2f2 x do. (3) 



where do, dP„, dp and SB are the perturbed velocity, to- 
tal pressure, density, and magnetic field, respectively. In 
terms of the perturbed gas pressure, one has 



dP„ = dP - 



B dB 

p 



(4) 



For our purposes, the linearized momentum equation is 
most useful in lagragian form. Introducing the lagragian 
displacement £, such that v = d£/dt, the continuity and 
induction equations can be integrated to yield 



*> = -V-(ft€), (5) 

dB = Vx(|xB ) 

= -B V-£-(£-V)B + (B -V)£, (6) 

while the expression of the total pressure perturbation fol- 
lows from our adiabatic equation of state: 



<5P, 



- € .VP - 7 J , V-€+^^ ) (7) 



SP 




where dP and dP m represent the gas and magnetic pres- 
sure perturbations, respectively. 

2.2. Ballooning formalism: 

The origin of MHD instabilities lies in inhomogeneous or 
non-static MHD equilibria. In this paper, we focus on in- 
stabilities related to equilibrium pressure-gradients, field 
line curvature and magnetic shear; they are amenable to 
an analytic description in the framework of the balloon- 
ing ordering (Newcomb, 1961). As this formalism is not 
well-known in the astrophysical community, we briefly de- 
scribe it below, and generalize it to rotating systems in a 
straightforward manner. 

The rationale of the ballooning expansion scheme fol- 
lows from properties of the linearized magnetic tension 



and of MHD wave propagation in homogeneous media. 
Let us consider plane waves, where £ oc exp(iu)t — ik ■ r), 
and assume that the direction of propagation is nearly 
perpendicular to the magnetic field, i.e. k /f <C fcj_ where 
k// and k± are the components of fc parallel and perpen- 
dicular to the magnetic field, respectively. In this limit, 
the frequencies of the slow [w| ~ CgV A /(Cg + V A )k 2 ] and 
Alfven = V^k 2 ] modes are very small compared to 
the fast one [uj 2 f ~ (C| + V|)fcjJ (see Appendix). 

On the other hand, inhomogeneous equilibria are char- 
acterized by a scale L a of inhomogeneity (e.g., \K C \ ~ 
\Kb\ ~ Lo 1 )- Order of magnitude considerations show 
that such inhomogeneities will contribute terms of order 
V^L~ 2 or CgL~ 2 to the dispersion relation, i.e. can desta- 
bilize the slow and Alfven modes inasmuch as the mode 
scale in the direction of the magnetic field, represented by 
kj x in homogeneous equilibria, is not significantly smaller 
than L (we leave aside the fast mode for the time being). 

With this consideration in mind, let us go back to wave 
propagation in homogeneous media, and introduce the or- 
thogonal reference frame (e /(/ , ej, &a) where e /7 = B j ' B Q is 
parallel to the unperturbed magnetic field, e; = kj_/k± is 
parallel to perpendicular component of fc, and &a = % x e; 
(see Fig. |l|); the subscripts I and A stand for longitudinal 
and Alfvenic, respectively (e /h e;, and are the direc- 
tions of the displacement of purely fast, slow and Alfvenic 
modes respectively in the limit of nearly transverse prop- 
agation adopted here). Denoting (£//,£;, £a) the compo- 
nents of the lagragian displacement £ in this reference 
frame, and introducing the small parameter e = k///k±, 
the momentum equation in the direction of e; implies 
that ~ for the slow mode; this implies in turn that 
<5P* = to leading order in e (see Appendix). This re- 
flects the fact that in quasi-perpendicular propagation, 
slow magnetosonic perturbations behave quasistatically 
and purely compressively in the fast direction (parallel 
to ei). This cancellation of the total pressure is essential 
from a technical point of view: it allows us to introduce a 
WKB-like approximation perpendicular to the magnetic 
field, because it can be used to remove all spatial deriva- 
tives of the mode amplitude in the direction perpendicu- 
lar to field lines, as argued in the Appendix; this greatly 
reduces the complexity of the problem. Incidentally, this 
eliminates fast modes from the analysis; in any case, these 
modes cannot be destabilized by the process considered 
here in this WKB-like approximation, as the terms con- 
nected to inhomogeneities are always subdominant in their 
dispersion relation. 

In this discussion, we have focused on instability, but 
the formalism can account for overstability as well, which 
can arise when considering non-static equilibria. It turns 
out that this possibility plays little role for the modes we 
consider (see next section); the question of overstability 
for the fast magnetosonic mode is considered in section 4. 

Turning back to inhomogeneous media, let us con- 
sider lagrangian displacements such that £ = £ Q exp(iujt — 
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field lines stretching 
«a fc- 




Fig. 1. Definition of the reference frame {e /h ei, ca) 



Fig. 2. Magnetic flux tube displacement 



i<£>). The standard ballooning formalism (e.g. Dewar and 
Glasser 1983) assumes k± = V$ _L B everywhere, intro- 
duces a small parameter e = (k±L Q )~ 1 , so that variations 
in the direction perpendicular to the unperturbed mag- 
netic field are treated in the WKB approximation, but 
not parallel ones. As suggested by the properties of wave 
propagation in homogeneous media, one introduces the 
following self-consistent ordering: 
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0{V\L-\ 



(8) 
(9) 



As shown in the Appendix, the longitudinal compo- 
nent of the momentum equation implies 5P* = to leading 
order in e, as expected, while the two remaining compo- 
nents reduce to: 



d 2 U 
dt 2 

dt 2 



'D a -U = Ca-£ j ,+Fa{£), 



(10) 



(11) 



where T stands for the inertial and Coriolis forces due 
to the rotation of the jet. In writing down these equa- 
tions, we have defined auxiliary quantities, C and T>, which 
characterize the perturbation of the magnetic tension; C 
introduces coupling between the Alfvenic and slow mag- 
netosonic part of the perturbation, and vanishes for ho- 
mogeneous media. Defining d/ds = e//.V, these quantities 
can be expressed as follows: 



V A -U = Vl J -^U + Kb, ^ 



-r. . +iCb„T. 



(rV) 



2 ^ CA 1 + (^ PA ~ ^ CA _ ^ bA 



Uh ( 12 ) 



'9£ 



for the Alfvenic equation, and 



Qt 



VI 



l__d 

fh ds 



Po Vg u (jC PA - JC CA - lC bA ^A 
(K b - K c )-T -£, A + (jC bA - !C CA ^j -J^- 



(13) 



(14) 



(15) 



for the slow magnetosonic equation; V S 2 M = 8Va/{1 + P) 
is the slow magnetosonic speed. The matrices T and T 
arc defined by 



= [(e // -V)e Q ± (e Q -V)e 7/ ] 



(16) 



a— {/,a} 



In homogeneous media, only the first term of T>a 
and V// remain. In this case, Eqs. ( ]l0| ) and ( pd| ) reduce 
to the dispersion relation of the Alfven wave and slow 
magnetosonic wave, respectively. The coupling coefficient 
C A cx IC CA , so that the coupling of the slow magnetosonic 
wave to the Alfven wave is directly related to field line 
curvature in the absence of rotation. The matrices l en- 
capsulate purely geometric effects due to the orientations 
of the field lines and of the wavevector k. 

These equations are usually written in a more compact 
form (see e.g. Dewar and Glasser 1983), which can be re- 
covered after some lengthy algebra by inserting the equi- 
librium relation for non-rotating jets and rescaling the dis- 
placement components as to include the purely geometric 
effects just pointed out. The expanded version presented 
here is more suitable for our purposes. 

Generically, pressure-driven instabilities work on a 
combination of the alfvenic and slow magnetosonic dis- 
placements, however, one can devise situations where one 
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or the other dominates. As an example, let us consider the 
origin of the destabilisation of the Alfven mode in a static 
shearless two-dimensional equilibrium with no variation 
along field lines (it is also possible to devise situations 
where one primarily destabilizes the magnetosonic mode, 
but the reasoning is less direct). Let us chose &a in the 
direction of the field line curvature, and consider displace- 
ments of a flux tube in the alfvenic direction, and constant 
along field lines [see fig. (0)]. In this approximation, only 
Eq. (|lo| ) remains, and reduces to 



dt 2 



4 = vU 2I31C CA K PA - aJ^K 2 



1 + 



(17) 



Instability follows if the right-hand side member of this 
equation amplifies the motion, i.e. if 



1C PA > 



This criterion reduces to the well-known criterion of 
instability with respect to m — modes (known as the 
"saussage" mode, which is a pressure-driven mode) in 
cylindrical systems, (Kadomtsev 1966; see next subsec- 
tion); it can be more directly established in the following 
way. 

Because the total pressure perturbation cancels, stabil- 
ity with respect to the applied displacement is controlled 
by the sign of the magnetic tension, which, in this case, 
reduces to 5Ta — 2p V^Sb // IC CA because the perpendicular 
perturbation of the magnetic field vanishes for such dis- 
placements. In this type of displacement, V.£ = V.£j_. 
The considered displacement of the flux tube is accompa- 
nied by a compression or an expansion (i.e. a change in 
V.£j_), whose sign and magnitude depends on properties 
of the considered equilibrium, and which can be computed 
from Eq. (A. 7), combined with Eqs. (||) and (fij). This pro- 
cedure finally yields 8b// = (i[K, PA — 2fC CA /{\ + (3)]£,a, whose 
sign leads to a destabilizing contribution of the magnetic 
tension when the criterion ( p"8| ) is satisfied. This reason- 
ing is typical of an interchange instability, and in fact, 
the ballooning formalism applies to both interchange and 
ballooning modes. 

2.3. Reduction to cylindrical symmetry 

The stability properties of the ballooning system of equa- 
tion has been been largely studied for cylindrical static 
equilibria in the context of thermonuclear fusion research. 
In particular, the well-known Suydam criterion provides a 
necessary condition for stability of the ballooning modes 
(see Freidberg 1987, chapter 10 and references therein). 

In this paper, we follow a different route by relaxing the 
eikonal condition fc^ = V$. Instead, we Fourier transform 
the displacement vector in the vertical and azimuthal di- 
rections: £ = £o(f) exp^iut—imO—ikz), i.e. we consider in- 
terchange, rather than ballooning modes. Ballooning and 



interchange stability analyses in static media based on the 
form of the energy principle derived by Furth et al. 1965 
show that the destabilizing term is cx K- p .K, c . For cylin- 
drical configurations, this is maximized by chosing in 
the radial direction; this choice also leads to the simplest 
calculations. Therefore we assume = — e r from now on 
(the minus sign is arbitrary, and chosen so that important 
quantities appearing later are positive). This completely 
specifies the reference frame (e /7 , ej, e^) introduced previ- 
ously. The existence of modes which self-consistently sat- 
isfy these constraints in the framework of the ballooning 
expansion is verified in section 3. 

Let us define k = meg/r + ke z and k// — k.B /B . 
The constraint of quasi-perpendicular propagation im- 
plies that k// <C k.ei; this is satisfied in particular in 
the vicinity of magnetic resonances, implicitely defined by 
mBg/r + kB z = (so that k /f = on resonant surfaces), 
where Bg and B z are the azimuthal and vertical compo- 
nents of the unperturbed magnetic field. In general, due 
to magnetic shear, the constraint of quasi-perpendicular 
propagation implies that our system of equations is only 
valid within a region of limited radial extent for any given 
mode. However, we are mostly interested in generic con- 
ditions for instability rather than on mode structure; this 
makes this restriction of little importance, since u> is con- 
stant throughout the radial extent of the mode, and not 
only in the region of validity of our system of equations, 
so that conditions over to established in this region are 
nonetheless valid for the whole mode. Furthermore, one 
can define a resonant k for any m and for an arbitrary 
magnetic surface, so that our analysis applies to the whole 
structure. Finally, Goedbloed and Sakanaka (1974) have 
proved an important theorem which implies that, at fixed 
m and k, if there is a (radially) local unstable mode, there 
must be a globally unstable one with maximal growth rate 
for this choice of m and k; this theorem applies to static 
equilibrium configurations, i.e. in our case, to nonrotating 
cylindrical columns. Its relevance for real jets is unclear, 
due to velocity shear and rotation, but it nevertheless sug- 
gests that local stability results (in radius) might be rele- 
vant outside their strict domain of validity. 

We conclude from this discussion that although the 
magnetic shear limits in theory the applicability of our 
analysis to a finite radial region for each mode, this limi- 
tation is of limited practical importance. Note furthermore 
that we have made no assumption concerning the radial 
behavior of the modes we consider, besides the limitation 
imposed by the nearly transverse propagation constraint, 
which in general should limit both amplitude and phase 
variations to scales much larger than kj 1 (however, see 
the discussion of section 3 for an important exception to 
this rule). 

With our choice of reference frame, and the help of 
the equilibrium relation Eq. (^), the system of equations 
derived in the preceding section reduces to 
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3.1. Nonrotating jets 



- i I 2V 2 



SM 



K CA k// 



2J2 



1+0 



+ 2Cl up ^-]L.\. (1!») 



[u; 2 - [k] + k 2 + k 2 r )V 2 \ U = 

2Vi M )C CA k // + T fp!f+2n o u^)$ / . (20) 
In these equations, 



In order to get some grasp of the meaning of Eqs. (|l9j) and 
(p0|), let us first consider non- rotating jets. In this case, the 
MHD perturbation equations become self-adjoint with our 
approximations, so that uj 2 is real. 

Interesting general constraints can be derived from the 
dispersion relation which reads 



^-[${vl + vl M ) + klvl]J 



k 2 (k 2 



kl)vl M vl = 0, 



(25) 



2 ^o r 

K = (i ZjZ °j3)v 2 CA + ^ bA ~ ^ PA ^ 



n 2 r 2 \ 1 



2K C . 



l + 0\V 2 



(1 + 0)V% r 2 



1 - 



(1 + 0)JC 



2K n 



(21) 



k a is defined in terms of the following auxiliary quantities 



where k 2 = TT + 20JC CA JC PA = 4/9X^/(1 + 0) - k 2 Q . The 
roots can readily be extracted from this dispersion rela- 
tion, but their properties are more directly understood in 
the following way. 

The product of the two roots is positive when k 2 > 
k 2 . One can check that this implies that the sum is also 
positive, so that both roots are positive. Therefore 



(26) 



and reads 
k 2 



r + r" 



23 



2IC 2 i 

CA \ B 2 



1+0 
B 2 \ fdluh 



h _ B z 
2jrr ~ B e ' 



2 B 2 B 2 rh' 
^~B± ft"' 



(JC CA +K,b A -K, PA )K CA . 



dim 



1 + 



1 - 



2K, C , 



2/3 K c 
1 + r 



(22) 



(23) 



(24) 



In Eq. (po|) , V^(fc 2 + k 2 ) plays the role of a general- 
ized Brunt- Vaisala, frequency; k r arises from the entrap- 
ment inertial force term. The term T T symbolizes the 
remaining component of the product of T matrices, and 
represents the effect of the magnetic shear (note that 
rh' '/h — dln/i/dlnr is the magnetic shear). The terms 
on the right hand side of Eqs. ( |l9| ) and ( |20| ) represent the 
coupling between the two modes due to field line curva- 
ture, the entrainment inertial force, and the Coriolis force, 
respectively. Note that JC CA = (Bg/B ) 2 jr. 



is a necessary condition of instability for the modes con- 
sidered here; it is also a sufficient condition, because the 
root product is negative then, implying that one root is 
negative. This relation implies in particular k 2 > 0, which, 
in the absence of magnetic shear, reduces to the necessary 
condition that the ballooning driving term in the energy 
principle of Furth et al. (1965) be destabilizing and is iden- 
tical to @. 

Another condition is obtained by taking the limit kg — > 
in the dispersion relation, which yields uj 2 — and uj 2 = 
as the two roots. This implies that 



v\kl 



k 2 <0 



(27) 



is also sufficient (but not necessary) condition of instabil- 
ity of the modes we consider. 

To pinpoint the meaning of the second criterion, it 
is instructive to examin the behavior of the modes when 
\k//\ <C \k c \, \k a \ (which can always be achieved by an 
appropriate choice of k for any particular value of m). 
The first mode has uj 2 ~ 0. Then, Eq. (|o|) yields 
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'2i0K, CA kf£j/ 



(l + 0)k 2 ' 
which can be inserted in Eq. (R~9|) to obtain 



(28) 



3. Dispersion relation and instability conditions 

The complete dispersion relation can easily be obtained 
from Eqs. ([l9]) and (po|). However, we are mostly interested 
in the fc /(/ — » limit, for reasons exposed below. 



UJ 



-V? 



SM' 



h 2 h 2 

w 



(29) 



For natural reasons, we refer to the modes which exhibit 
this behavior in the fe» — > limit to slow (magnetosonic) 
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modes. From the general solution to the dispersion rela- 
tion, the maximum value of uj 2 for this mode is obtained 
for kj ~ k 2 , which provides us with an estimate of their 



maximum growth rate, once inserted in 
mode has 



vik 2 . 



The second 



(30) 



Then Eq. <M) implies 



(31) 



-2z/3/C CA fc//^ 

{i+m ■ 

We refer to modes which exhibit this behavior in the k f/ — > 
limit as to Alfven modes. The maximum growth rate for 
these modes is reached for k f/ =0. 

With these definitions, the slow modes are unstable 
and the Alfven ones stable when k 2 > k 2 and k 2 > 0; 
the situation is reversed when k 2 > k 2 and k 2 , < 0. Both 
modes are stabilized by the usual contribution to the mag- 
netic tension when k 2 > k 2 . 

Note that when the magnetic shear is important, £ for 
the Alfven mode, and £4 for the slow mode, undergo fast 
radial variations. One can check that fast radial variations 
of £j/ are not inconsistent with our procedure, which is not 
true for £4. Therefore, when this occurs, one must discard 
the slow mode, which is not a self-consistent solution of 
our system of equations. 



3.2. Jets with vanishing current density on the axis 

There are situations in which one expects that Bg in- 
creases with r faster than r in jet inner regions (in 
which case the current density vanishes on the axis). In- 
deed, for jets which are launched from an accretion disk, 
u z Bg/B z — r(f2 — f2+) where fi and are the matter 
and field rotation respectively, u z is the jet velocity along 
the jet axis (Pelletier and Pudritz, 1992). In these regions, 
one usually expects u z ^ 0, while Q, 57* — > when r — ► 0. 
Furthermore, polarisation observations indicate that the 
magnetic field is predominantly aligned with the jet axis 
in jet cores, at least for quasar jets (see Gabuzda, 1997, 
and references therein). It is therefore reasonable to ex- 
pect that \Bg\ < or <C \B Z \ over a sizeable fraction of 
the inner region of these jets. In these regions, fcg < 
is automatically satisfied because the shear term domi- 
nates (as rotation is unimportant there) and is negative 
[see Eq. (|2"i|)], so that the inner regions of these jets are 
generically unstable. Note that our neglect of the velocity 
shear in the derivation of our equations has no impact on 
this conclusion, because both rotation and velocity shear 
terms are expected to be negligible in the considered jet 
regions. Obviously, \B Z \ must exceed \Bg\ over a sizeable 
fraction of the jet interior for the instability to have no- 
ticeable effects and growth rates. 



The destabilizing action of the magnetic shear just 
pointed out is the most important finding of this paper. 
It can be given a heuristic explanation in the following 
way. Regions with a destabilizing contribution of the mag- 
netic shear are regions in which the shear contributes to 
an extra-increase of the Alfvenic component of the equilib- 
rium tension when one moves outwards. Applying a purely 
radial displacement to a magnetic flux tube (i.e., consid- 
ering an Alfvenic mode with k ff = 0) produces a varia- 
tion of the Alfvenic component of the magnetic tension 
which does not include this extra piece due to the mag- 
netic shear; therefore, the variation in the magnetic ten- 
sion produces a restoring force on the displaced flux tube 
which is smaller than the restoring force on the equilib- 
rium flux tube at the same location, and instability fol- 
lows. 

Note also that we must discard the magnetosonic 
mode, which is no longer a valid solution of our equa- 
tions due to the fast amplitude gradients of £4 induced by 
the shear, as explained in the preceding section. 

3.3. Nonrotating jets dominated by the toroidal field 

In magnetically confined jets, one usually assumes \Bg\ 3> 
\B Z \ in the confining region ("Z-pinch" configurations). 
This follows in particular when jets are launched from ac- 
cretion disks, because they must open considerably before 
the critical surfaces are crossed, and the magnetic colli- 
mation becomes important. This boosts the ratio \Bg/B z \ 
first because the poloidal flux is conserved within magnetic 
surfaces, second because the toroidal field is considerably 
stretched in the process, while this ratio is expected to be 
of order unity on the disk "surface" (Ferreira 1997; Casse 
and Ferreira 2000). In this case, k// ~ m/r + kB z j\Bg\ and 
k± ~ k. Furthermore k// depends on the radial location; 
Eq. ( |3l| ) shows that £ contain no large amplitude varia- 
tions [as implicitly assumed in the derivation of Eqs. (19) 
and @] provided that \krB z /B e \ < 1 and |m| < 1 (note 
however that fast amplitude variations of £a do not inval- 
idate the derivation of these equations). Then, the con- 
dition of quasi-perpendicular propagation is satisfied for 
I Air I 3> 1 (as m is a small integer). 

When these conditions of consistency of the slow mode 
with our approximations hold, one can eliminate all gradi- 
ents in terms of logarithmic gradients of Bg (the gradient 
of B z vanishes or is negligible). This allows us to reexpress 
k as 



1 + /3 



1 



20 



f3\ ( d\nrB 6 
r- 



dr 



(32) 



The shear term has disappeared because it is negligi- 
ble for I I ^> \B Z \; the first term in brackets represents 
the stabilizing action of the plasma compression (it arises 
from the divergence of the displacement), while the other 
term represents the action of pressure gradients, either 
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stabilizing or destabilizing. The expression for k 2 is very 
similar, and obtained by changing the sign and removing 
the compression term. 

The limit \Bg\ > \B Z \, completed with \kr\ > 1 and 
\krBg/B z \ < 1 allows us to connect with previously de- 
rived results on pressure-driven instabilities. In this case, 
one has k f/ = sgn(Bg)m/r + kB z /\Bg\, JC CA = 1/r, and the 
criterion (5ft) becomes, with the help of Eq. (|35~ 



d\n\Bg\ 1 



fcr^l -1. 



(33) 



which has been derived for m = 1, k = modes with 
V.£ = by Kadomtsev (1966), and for an arbitrary mode 
by Begelman (1998). 

The condition ( pTj ) which identifies the Alfven branch 
as the unstable one reads in the same limit: 



It is useful to characterize the jet rotation by intro- 
ducing x = ^o 7 " 2 V^a- The equilibrium condition Eq. ^ 
implies that x, like /3, is of order unity at most. More 
specifically, one has (Pelletier and Pudritz 1992) 



X 



1 (1 



M 2 r\/r 2 ) 2 



M 2 (l-r 2 A /r 2 ) 



(35) 



where M is the poloidal Mach number and is the Alfven 
radius. In strong jets, one has M 2 r 2 A /r 2 ~ 1 whereas r is 
always larger (possibly much larger) than ta, which, as 
stated, implies that x is of order unity at most. We assume 
X ;$ 1, which maximizes the influence of rotation on the 
dynamics. 

The influence of the jet rotation is again analyzed by 
deriving the dispersion relation, which reads: 



d\n\Bg\ 8 
dlnr > J- 



r2 , y2 



(34) 



a condition also obtained by Kadomtsev (1966) and Begel- 
man (1998) for the m — "saussage" mode (remember 
that our definition of B differs from the standard one by 
a factor 7/2). 

We point out again that (| 



|) [or equivalently (J26J)] is 
a necessary and sufficient condition for having unstable 
modes, but (34) is a condition for the Alfven branch to 
be unstable; these conditions do not assume any partic- 
ular value of m and do not rely on an incompressibility 
assumption. These points are often overlooked in the lit- 
erature. Note however that for smooth profiles, k 2 > 1/r 2 , 
so that ([33]) limits m to small values, independently of the 
limit obtained from the consistency requirement of the 
slow mode. It is well-known that condition (|34|) is easy to 
meet in practice; in particular, it is always satisfied for a 
Lorentzian (or Bennett) profile. ) We conclude with Begel- 
man (1998) that the configurations usually considered for 
magnetically confined jets are generically unstable in the 
absence of rotation. Indeed, the condition for instability 
Eq. (|2^) is always satisfied for m = 0, and probably m = 1 
as well in jet inner regions, because \Bg\ must increase 
from zero on the axis to values \Bg\ ^> \B Z \ in inner re- 
gions; in general it decreases too fast in most jet models 
when moving towards the jet radial boundary for outer 
regions to be unstable. Furthermore Eq. ( |34j ) implies that 
the Alfven mode is very likely to be the unstable mode. 



SMJ 



-4n o {B /B o )k // Vg M [2IC 



(k 2 + k 2 M + m 2 (B 2 /B 2 )] co 2 
X 



Br 



UJ 



k» k c — - k T 



(1 + 0) >2/V " 



X_ 

0r 



VImVI = 0, 



(36) 



where k 2 and k 2 are given in Eqs. ( pi]) and (|2l|). 

In the limit k// = 0, this equation has two roots, the 
magnetosonic branch (whose degeneracy is broken by ro- 
tation) uj 2 = 0, and the Alfvenic branch 



uj 2 = Vl 



l<>- I ^± 

'■' B 2 



(37) 



The last term in this relation is due to the Coriolis force 
and is always stabilizing. This is a natural consequence of 
the form of this force, f2 x £, which gives rise to coupled 
oscillations between the components of the displacement 
vector. The factor Bg/B Q arises because of the inclination 
of the magnetic vector B Q with respect to the rotation 
axis. 

The discussion of the stability properties in the ab- 
sence of rotation implies that stabilization by rotation will 
only occur when % ~ 1. To quantify somewhat this state- 
ment, let us rewrite k 2 and k 2 as 



3.4- Rotating jets dominated by the toroidal field 

Rotation, in particular the Coriolis force, introduces an 
extra-coupling between the Alfvenic and magnetosonic 
part of the displacements, which is not vanishing in the 
k/i — > limit. Nevertheless, we stick with our distinction 
between Alfven and slow magnetosonic modes, although 
it is less meaningful in this context. 



ki = 



1 4/3 

V 2 1 + B 



U2 _ 2 X 
r r 2 (l + B) 



1 + Bd\nr\Bg 
~2/3 dlnr 



l + Bd\nr\Bg 



2B d\nr 



2X 



i a+py 



(38) 



r 2 B(l + B)' 



(39) 



Evy Kersale et al.: Pressure- and magnetic shear- driven instabilities in rotating MHD jets 



9 



It is apparent that the contribution of the jet rotation 
to k 2 is always stabilizing. Furthermore, the pressure con- 
tributions to k 2 , and k 2 are destabilizing when condition 
( |34| ) is satisfied. The arguments of the previous subsection 
make this very likely to happen in the inner region of jets. 
A condition on the critical gradient of Bg can easily be 
extracted, but it is not very informative. Instead, limits 
of instability can be obtained in the following way. We 
assume that d\n\Bg\/d\nr ~ 1; this applies in jet inner 
regions when the behavior of the magnetic shear is not 
controlled by rotation, as in the previous section. The dis- 
persion relation (|3^) implies then that (3 < 1/4 for x — 1 
or reversely x < 1/2 for /? ~ 1 constitute sufficient con- 
ditions of instability. This implies that it is very difficult 
in practice to stabilize the Alfven mode by rotation, espe- 
cially that jets possessing an important cylindrical asymp- 
totic limit have \ ^ 1 from Eq. (|3^), as such jets have 
M > 1. 

Note that stability occurs for high (3, whereas in the 
fusion context, the reverse is true. This follows because 
destabilisation is provided by the rotation entrainment 
force rather than pressure forces. 

The stability properties of the slow mode can be dis- 
cussed in the k f/ — ► limit. Assuming that the frequency 
of magnetosonic modes u> oc k// for small enough k /h as 
implied by the dispersion relation, and keeping only the 
leading order term in k /; yields a reduced dispersion rela- 
tion: 

Alu 2 + Blu + C = 0, (40) 

where 

A = [(k 2 Q +k 2 )V 2 + AQ 2 (B 2 /B 2 )], 
B = m {B g /B )k ff Vl M [2K CA + 

c = t (-*2 + ^ - ^ ( 2/c - + £) ) v ^ 

(41) 

Because of the term oc lo in this reduced dispersion 
relation, overstability rather than instability occurs when 
the discriminant A = B 2 — A AC < 0. It turns out that for 
X — 1 and d In Bg jd In r ~ 1 , A > whatever the value of 
/3, implying that the slow mode is stabilized by rotation 
in the small k /f limit. 

This suggests that the slow mode is more easily sta- 
bilized than the Alfvenic one. In any case, our result on 
the Alfven mode is sufficient to conclude that jets with a 
dominant toroidal field are difficult to stabilize by rota- 
tion. 



4. Discussion 

In this paper, we have focused on the stability of rotating 
MHD jets with respect to inhomogeneities connected to 
pressure gradients and magnetic shear, in the framework 
of the ballooning ordering. The most remarkable finding 
of our analysis is the possible destabilizing influence of 
the magnetic shear pointed out in section 3. Indeed, the 
shear appears quadratically in Suydam's well-known crite- 
rion, which applies to the stability of cylindrical columns, 
whereas the stability criterion of the preceding section is 
linear in magnetic shear. We point out that, had we not 
relaxed the eikonal condition in section 2, we would have 
recovered Suydam's criterion (see Freidberg 1987, Chap- 
ter 10). This unexpected behavior of the magnetic shear 
implies that the inner regions of jets launched from ac- 
cretion disks are generically unstable, as argued in section 
3. 

We have also recovered that for jet models dominated 
by a toroidal magnetic field, jet interiors are likely to be 
subject to this type of instability, a result already pointed 
out by Begelman (1998). Furthermore, we find that rota- 
tion has little influence on this conclusion; rotation can 
only stabilize jets when it is close to the limit allowed by 
the radial equilibrium condition, and for values of (3 close 
to the same limit. 

Our results on pressure-driven instabilities agree with 
those found in the literature, most notably Begelman 
(1998) and Ferriere et al. (1999), once appropriate limits 
are taken. Begelman's study is restricted to Bg/B z 3> 1, 
which removes the effect of the magnetic shear from the 
analysis. Ferriere et. al. (1999) do not rely on the bal- 
looning expansion, keep the full six-dimensional disper- 
sion relation, and discuss in detail the effects of rota- 
tion and gravity, besides gas pressure and field curvarture. 
They give a detailed and interesting synthesis of gravita- 
tional and centrifugal instabilities in magnetospheric stud- 
ies. However, to achieve such a generality, they rely on a 
particular coordinate representation of the magnetic field, 
which is well-adapted to planetary magnetospheres, but 
forbids also the presence of magnetic shear. 

Jet models which are not dominated by the toroidal 
magnetic field have been proposed in the literature. One 
example is the force-free jet model of Appl and Camen- 
zind, 1992, which has Bg ~ B z , and belongs to a class of 
configurations known as "reversed field pinch" in the fu- 
sion research context. These configurations are known to 
be stable with respect to pressure-driven instabilities for 
moderately low values of j3 (and, more generally, to pos- 
sess good stability properties). This apparently conflicts 
with our results on destabilization by the shear. However, 
we point out that destabilization occurs because Bg varies 
faster than r in the jet inner region. This variation is in- 
duced by the behavior of the jet matter and field rotation, 
a feature which is absent from fusion control devices. 



10 



Evy Kersale et al.: Pressure- and magnetic shear- driven instabilities in rotating MHD jets 



In this analysis as in all others we are aware of, we have 
assumed that the jet equilibrium structure is invariant 
along the jet axis. This need not be the case, as exampli- 
fied in various self-similar solutions (Sauty and Tsinganos, 
1994; Contopoulos and Lovelace, 1994), and more recenly 
in a complete analytic solution displaying a more com- 
plex magnetic field structure (Bogoyalenskij , 2000). The 
relevance of the available stability analyses to such con- 
figurations is unclear; furthermore, these configurations 
are possibly unstable with respect to appropriately cho- 
sen ballooning modes. In any case, this question should 
be addressed before firm conclusions can be drawn. 

We have not tried to derive precise growth rates, as 
these depend on the radial structure of the most unstable 
modes, the determination of which is outside the scope 
of the approximation scheme we have adopted. However, 
order of magnitude considerations show that these insta- 
bilities should develop within a few dynamical time-scale 
(~ aV^ 1 where a is the characteristic radial size of the 
considered region). 

Let us also return to the question of overstability, 
which was briefly alluded to in section 2 and 3. Our dis- 
persion relation shows that lo 2 is real for alfvenic modes 
in the k a — limit, implying that lo is real when lo 2 > 0. 
This property remains true when to(k//) is computed as a 
scries in k ff , to all orders in this series if it is true to low- 
est order (this can be shown by induction from the full 
dispersion relation, because it involves only real terms). 
Therefore the stability property of the Alfven mode is ob- 
tained in the k f/ = limit. In the same way, the dispersion 
relation of Ferriere et. al. (1999) [their Eq. (27)] implies 
that the fast magnetosonic mode is stable in the absence 
of magnetic shear in the perpendicular WKB approxima- 
tion. This suggests that the fast mode is probably not 
destabilized in the asymptotic limit examined in this pa- 
per. 

We point out again that, in spite of their partly local 
nature, pressure-driven instabilities, and magnetic shear- 
driven ones as well, are likely to produce large-scale dis- 
ruptions of MHD jets. How real jets survive these instabil- 
ities as well as Kelvin- Hclmholtz and current-driven insta- 
bilities is a still open and pressing problem, especially that 
3D instabilities seem more vigorous than there 2D coun- 
terparts (see, e.g., Bodo et al, 1998). However, progress on 
this topic has recently been accomplished. For the Kelvin- 
Hclmholtz instability, recent 3D simulations have shown 
that for initial magnetic Reynolds numbers smaller than 
~ 50, the shear layer broadens, the turbulence level de- 
creases and the flow becomes quasi-laminar and stable 
in the end (Ryu et al., 2000); as jets are not expected 
to be drastically out of equipartition, and as the Lorentz 
force provides the origin of their acceleration, this proba- 
bly ensures their stability against a possible disruption by 
the Kelvin-Hclmholtz instability. Also, the kink instabil- 
ity seems to lead to helicoidal equilibria with a redistribu- 
tion of the current density rather than to disruption (Lery 



et al., 2000). Finally, it has recently been recognized in 
the fusion literature that pressure-driven instabilities can 
be stabilized by large equilibrium gradients, a somewhat 
counterintuitive result (Terry, 2000). 

One of the interesting questions raised by this work 
concerns the possible link between the type of instabili- 
ties we consider and the particle acceleration mechanism 
which is responsible for the high energy emission seen 
in AGN jet interiors. For example, in the context of the 
two-flow model, most of the high energy emission comes 
from a highly relativistic pair plasma occupying the in- 
ner region of a slower jet which contains most of the 
mass and provides most of the energy. This model re- 
lies on energy injection to the pair plasma from the bulk 
of the jet through turbulence (Renaud and Henri, 1998); 
note also that the existence of a pair plasma in AGN jets 
has recently received an interesting observational support 
(Wardle, et. al, 1998; Hirotani et. al, 1998). In such a 
picture of the high energy emission, the required turbu- 
lence would most likely be provided by MHD instabilities 
in the nonlinear regime, as other sources look unpromis- 
ing. Indeed, a purely hydrodynamical origin in regions of 
increasing angular velocity with radius, as advocated by 
Richard and Zahn (1999) is not found in numerical sim- 
ulations (Balbus et al., 1996), and happens in laboratory 
Couette experiments (performed in general with narrow 
gap widths between the inner and outer cylinder) only 
when the outer cylinder completely dominates the rota- 
tion (Coles, 1965), a situation which has little connection 
with disks and jets; furthermore, the magneto-rotational 
instability (Balbus and Hawley, 1991) does not work in 
this regime, and the velocity shear is ineffective in driving 
the Kelvin- Hclmholtz instability in the subalfvenic portion 
of the jet in jet interiors. 

To conclude, let us finally point out that pressure- 
driven modes are also known to be avalanche-like instabil- 
ities, due to the existence of a critical gradient, and this 
feature raises important question concerning the action of 
this instability on jet profiles, as well as their role in jet 
variability. 

Appendix A: MHD perturbation equations and 
the ballooning approximation 

Let us first recall some properties of wave propagation in 
homogeneous media. In the reference frame introduced in 
section 2.2, the momentum equation yields the following 
three component equations 

{u 2 -C 2 k 2 ) ^ = C 2 s k„k ± ^, (A.l) 
(lo 2 -C 2 k\- V% k 2 ) 6 = C 2 k„ k ± , (A.2) 
(u; 2 -V 2 k 2 ) U = 0, (A.3) 
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while the total pressure perturbation becomes 

dP, = -i Po [(C 2 S + V%) fcx6 + C% k,(,] . (A.4) 

Eq. ( |A.3| ) gives the dispersion relation of Alfven waves, 
lo\ — V\kf n which decouple from the two magnetosonic 
modes described by the remaining two equations. In the 
approximation of interest here (e = k^/k± <C 1, i.e. 
nearly perpendicular propagation), these two equations 
imply u\ ~ C 2 s Vl/{Cl + V%)h* and & ~ 0(e&) for the 
slow magnetosonic wave, while up ~ (C| + V^)k« and 
£// ~ 0(e£ i i) for the fast magnetosonic one. 



Furthermore, the momentum component Eq. (A. 2) 
combined with Eq. ( |A.3| ) implies that 5P* = to leading 
order in e. 

Let us now turn to inhomogeneous and non-static me- 
dia. The ordering assumed in Eqs. (||) and implies that 
the longitudinal component of the momentum equation 
[Eq. (||)] reduces to 



6P* = SP 



B -5B 



= 0. 



(A.5) 



to leading order in e. In establishing this result, the un- 
perturbed momentum equation [Eq. (|^)] has been used to 
show that Q ~ u at most. The total pressure perturba- 
tion cancels because in the quasi-perpendicular propaga- 
tion assumed here, slow and alfvenic perturbations behave 
quasistatically in the fast direction, and because fast per- 
turbations are only compressional in the fast direction in 
homogeneous media, implying that inhomogeneous cor- 
rections are negligible in our WKB-likc approximations. 
To derive the remaining two momentum component equa- 
tions, we need to evaluate the tension perturbation in 
terms of the displacement £ [the entrainment inertial term 
follows immediately from the perturbed density Eq. (|^)] 



ST = 



B VSB 5B-\7B Q 



/' 



(A.6) 



The magnetic field perturbation 8b = SB/B a [cf 
Eq. (Jsj) ] along and perpendicular to the unperturbed mag- 
netic field read 



A; 



-V-^-(7C hi +7C c )-a, 



In Eq. (A. 7), the first term represents the variation of 
the field in its original direction due to the conservation of 
magnetic flux in transverse compression of the plasma; the 
second term represents the contribution of field advection, 
the last one the contribution of field line stretching. 

With these results, the perpendicular component of 
the tension perturbation reads 



{l,A} 



PoVl 



— < 

ds 



>{l,A} 



{l,A} 



+ e {l , Ay T ■Sb x + 2Sb,K e{liA} ) , (A.9) 
while the parallel tension perturbation is given by 



sr // = PoV J 



d8j„ 



(A.10) 



These equations for the perturbed magnetic field and mag- 
netic tension are exact (the ordering has not been used in 
their derivation). Furthermore, by noting that 



ds 



6 IC b// - V-£, 



(A.ll) 



it appears that derivatives of the displacement in the mag- 
netic force term enter only through derivative along the 
field line and V.£; the same crucial feature holds for the 
entrainment inertial force term, while the derivatives of 
the displacement do not enter the expression of the Cori- 
olis force term. With the help of the cancellation of the 



total pressure perturbation Eq. (A.5), one can express V.£ 
in terms of £ and of the derivatives of the components of 
£ in the field line direction. This yields 



1 



1 + 13 

Ah 



(A.12) 



{1 + /3) (JC bA + ]C CA )U} 
This reduces the problem to the set of ordinary differential 



(A. 7) equations for ^ and £4 given in Eqs. (10) and (O 



»± = r -Cl + tt- e i + e A + {£-K. c ) e„ 



ds 



ds 



r„ 6 + r 1A u + 

r Al 6 + r AA ^A + 



ds 
dU 

ds I 



(A. 



where the matrices T and T are defined in Eq. (|l6|), 
d/ds = e^-V, and the various ICs are defined right be- 
fore Eq. @. 
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